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Abstract. The problem of the existence and uniqueness of the motion of 
the system of an arbitrary number linear oscillators under a generalized dry- 
friction type control is studied. This type of control arises in the problem of 
steering the system to equilibrium. The problem of existence and uniqueness 
of motion under the suggested control is resolved within the framework of the 
DiPernaLions theory of singular ordinary differential equations. 
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1. Introduction 

The content of the paper is based on the author’s talk at XII All-Russian con¬ 
ference on control problems [1]. B.T. Polyak have expressed a great interest to the 
author’s presentation. This is why we are especially pleased to publish this paper 
in the jubilee issue of the “Automation and Remote Control”. 

The control of oscillatory systems is one of the most important issues of the 
optimal control theory [2]. A classical achievement in this area is the explicit con¬ 
struction of the feedback control for the minimum-time damping of a single linear 
oscillator [3]. From the point of view of the canonical system of the maximum 
principle, the problem is completely integrable. 

The subject of the present paper is a more complicated problem of control for a 
system of an arbitrary number N of linear oscillators with different eigenfrequencies 

/^\ Xi = yi 

^ yi = -i^iXi+u, |m| < 1, i = l,...,N, 

which is probably not completely integrable. 

Authors do not expect to design an optimal control analytically, however the 
problem of damping of oscillations of system (1) is to be solved. A standard way to 
damp oscillation goes via friction. The straightforward application of the Coulomb 
law leads to the vector control Ui = — signi/i. The scalar control we are using is a 
generalization of the dry friction in the sense that it has the form 

N 

(2) u= - sign'^Xiyi, 

i=l 

where Xi are some coefficients. This type of control arises naturally in author’s 
papers [4-6] based on the study of asymptotic behavior of reachable sets of the 
system (1). 

It is well known that the use of dry friction, although it helps damping oscil¬ 
lations, does not necessarily leads to a complete stop of the system: there arise 
standstill zones where the system is not moving at all, although the terminal set 
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is not yet reached. In [4-6] an approach to design of quasi-optimal control based 
on a combination of three different strategies has been suggested. At low energies, 
the problem is resolved by using the common Lyapunov functions [6], whereas at 
high and intermediate energies the control used takes the form of the generalized 
dry friction. 

A characteristic feature of control problems is the discontinuity of the right-hand 
sides of the differential equations of motion. One encounters the same phenomenon 
in the study of the motion under the dry friction. The question of existence of the 
motion is traditionally resolved by using the Filippov theory of existence of solutions 
of differential inclusions [7]. However, the problem of uniqueness of solutions goes 
beyond the Filippov theory, although the intuitive concept of motion stipulates an 
unambiguous determination of the trajectories of the system by the control law. 

In this paper, the problem of existence and uniqueness of the motion of the 
system of oscillators under the action of a control in the dry friction from is in¬ 
vestigated. We explain how the control in the dry friction form arises from the 
structure of asymptotic behavior of reachable sets. In the present work, the central 
issue is the existence of uniquely defined motion under the control. This problem is 
resolved in the framework of the DiPerna-Lions theory of singular ordinary differ¬ 
ential equations (ODE) [8,9]. In this paper, we confine ourselves to the statement 
of the existence and uniqueness theorems. Proofs can be found in [4]. 

Consider a system of an arbitrary number of linear oscillators under a common 
bounded scalar control. The systems has the from: 

(3) X = Ax + Bu, X u s U = R, ju] < 1, 

where the matrix A and the vector B are as follows: 

(4) ly A = diag(A,), B = (BB,. 


It follows from the Bellman dynamic programming principle that the minimum¬ 
time control has the from [10]: 


( 5 ) 


i{x) = -sign{B,p{x)), 


dT 

p = p{x) = 


Here T{x) the optimal damping time, p is the outer normal to the reachable set 
from zero 'D{T{x)), whose boundary passes through x, and the angle brackets stand 
for standard scalar product in Recall that the reachable set V{T) is defined as 
the set of point reachable from zero at the time T along the trajectories of system 
(3). 

In the real life, the reachable sets are quite complicated. This is why we will design 
a control by using Eq. (5), where p is the normal to an approximate reachable set 
such that its boundary passes through x. According to the asymptotic theory of 
reachable sets, for system (3) a good approximation to ViT) as T —>■ oo is provided 
by a set of the form TO, where H is a time-independent convex body [11-13]. As 
it is well known, a closed convex body M is uniquely determined by its support 
function [14], 


( 6 ) 


Hm{C) = sup(^,x). 

x£M 


The language of support functions [11] is convenient for a precise statement on 
the asymptotic behavior of reachable sets: 


Theorem 1. Suppose that the momentum p is written in the form p = (pi), where 
Pi = if.i,Vi)> * = 1) ■ • ■; the dual variable for Xi, rji is the dual variable for 
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Ui, and Zi = {rjf+oj^ absence of resonances, i.e., non-trivial relations 

between eigenfrequencies of from 

N 

(7) ''^miUJi = 0, 0 ^ TO = (toi, ..., mAr) G 

1=1 

the support function Ht of the reachable set 'D{T) has as T —> oo the asymptotic 
form 


r2Tr 


r. 2 ir 


( 8 ) Ht{p) = 


(fZn) 


N 


N 


Zt cos (fit 


difi ... d(pN + o{T) = TSj{z)+o{T), 


and the support function of the closed compact Q is given by the main term Sj{z). 


1.1. Formula for the support function. The support function Hq{p) of the 
convex body is given by the main term of asymptotic equality (8): 

N 

ZiCOS(pi dip, z = {zi ,..., zn) e 

If TV = 1, we obtain Sj{z) = ^\z\. In the case of two oscillators N = 2, the 
function 


(9) 


Haip) = Sj{z) = J 


( 10 ) 


Sj{z) = J \zi COSipi+Z2 cos ip2\ dip 


can be expressed via elliptic integrals [4]. In general, function 


( 11 ) 


^^{z) = 


1 


(27r) 


N 


N 




2=1 


N 


\{{l-t1)-^/^dh...dtN 


i=l 


is an Euler-type integral that defines a generalized hypergeometric function in the 
sense of l.M. Gelfand. 


1.2. Equation defining the control. The vector p is normal to the boundary dO, 
at the point {p). Therefore, the normal vector to the approximate reachable set 
pVt such that its boundary passes through x can be defined via the equation 


( 12 ) 


1 dHn{p) dSj{z)dz 
dp ~ dz dp' 


where p G and p > 0 are unknown. We note that the support function Hq 
is differentiable, and equation (12) has only one solution due to the smoothness of 
the boundary of the set fl established in [12]. 


1.2.1. Duality. More generally, the relation between the support function H = 
Hq{p) and the function p = p{x) in equation (12) is similar to the Legendre trans¬ 
form: 


(13) {x,p) = p{x)H{p), p{x) = max (a:,p), H{p) = max {x,p), 

H{p)<l p{x)<l 


and the corresponding point transformation x ^ p has the form 
(14) X = p{x)^{p), p = H{p)^{x). 

Here p and x are maximum points in equations (13). These relations make sense 
if the functions H and p are norms of their arguments. The differentiation in (14) 
refers to taking subgradients. If one of the function H or p is strictly convex and 
differentiable outside zero, then the other is so. Then relations (14) hold in the 
classical sense. 
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By using the function p = p(x), one can rewrite control (5) in the form 

(15) u(a;) =-sign{^B, 

Equation (15) is used even in the resonant case, when the asymptotic equality (8) 
does not work. Note that u{x) is a multivalued function, because sign(O) can take 
any value in the interval [—1,1]. 


1.2.2. Hamiltonian structure. For the considered control, the function p = p{x) 
from equation (12) plays the same role as the optimal damping time T{x) for the 
optimal control. The “canonical” momentum in (12) is p = —dp/dx. The function 
p{x) is the norm of the vector x in the metric such that the body is the unit ball 
wrt to it. This is a smooth function outside zero. 

Furthermore, a “maximum principle” holds: the compound vector (a;,p), where 
p = dp/dx, satisfies the Hamiltonian system 


(16) X = AxBu, u = sign{B,p), p =-A*p+^^B sign(^B,^ 


where the Hamiltonian has the form 

(17) n = {Ax,p) + \{B,p)\- 


dx 


The Hamiltonian vanishes along trajectories of the controlled motion in complete 
similarity with the case of the optimal control. 


2. Motion under the control 


Formally speaking, the motion of the oscillator system under control (15) is 
described by the differential equation 


(18) 


X = Ax — B sign 



In fact, as it is already mentioned, the right-hand side of this equation is not defined 
uniquely (sign(O) = [—1,1]). Therefore, we are dealing with a differential inclusion. 

Fillippov theory [7] says that the Cauchy problem for differential inclusion (18) 
is solvable for any initial condition a:(0), i.e. there is a function of time x(t), which 
is absolutely continuous and satisfies (18) at points of differentiability. This follows 
from the fact that the right-hand side f{x) is, first, linearly bounded |/(a;)| < C'|a;|. 
Second, it takes convex and compact values. Finally, the function f{x) is continuous 
as a multivalued function: if G f{xn) and Xn —?• x, the y G f{x), where y is any 
limit point of the sequence ?/„. Generally speaking, the Filippov “motion” is not 
defined uniquely: several trajectories from a single starting point are possible. 

The central result of the paper is that the motion under control (15) can be 
defined uniquely. Toward this end, the DiPerna-Lions theory of singular ODE [8] 
is used. The general idea is to build a global phase flow, perhaps, not everywhere 
defined, instead of finding an individual solution of the Cauchy problem for each 
initial condition. 


2.1. DiPerna—Lions theory. If b{x) is a Lipschitz function, then the Cauchy 
problem for ODE 


(19) X = b{x), x(0) = xo 

and for the partial differential equation 


( 20 ) 


dv dv 

»=L'>-Wa;r. 


v{x, 0) = Vo(x) 
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are equivalent. Using the method of characteristics, one can show that the solution 
V of (20) is given by the formula 

(21) v{x,t) = vo{(l)t{x)), 
where (pt is the phase flow for (19). 

In the classical work of DiPerna and Lions [8], the Lipschitz condition on b is 
substantially relaxed, and it is shown that the solution of problem (20) still exists 
and is unique and is given by formula (21), where (pt is the measurable flow. Roughly 
speaking, it is shown that it suffices that the Lipschitz condition holds in the integral 
sense, not pointwise. 

Definition 1. A weak bounded solution v of the Cauchy problem is called renor¬ 
malized solution, if for any smooth function /3 : R ^ R function pi{v) is again a 
weak solution. 

Using this definition, we formulate the following theorem: 

Theorem 2. Suppose that (extended) DiPerna-Lions conditions hold in the Cauchy 
problem for the transport equation (20): 

(22) div6 GL-, 6 G = Wioc, + 

where BV\oc = is the space of functions such that their derivatives are locally 

finite measures. Then there exists a unique measurable flow (pt{x) such that if vq{x) 
is a bounded measurable function, then the function 

(23) v{x,f) = v{pt{x)) 

is the unique renormalized solution of the Cauchy problem (20). 

We note that the Lipschitz condition can be restated as db/dx G L°°. This 
theorem can be easily proven by methods of [9], where a strengthening of the 
results of [8] is suggested. 

Corollary 1. The Cauchy problem for the transport equation that corresponds 
to Hamiltonian system (16) and a bounded initial condition vo{x,p) has a unique 
renormalized solution v. The solution has the form 


(24) 

v{x,p,t) = vo{pt{x,p)), 

where 


(25) 

(Pt ■■ R'^^, t G R 


is a uniquely defined measurable flow. Each curve 

(26) t ^ {x{t),p{t)) = pt{x,p) 
is absolutely continuous and satisfies (16). 

Theorem 2 and its corollary are useful, but they do not define any flow in the 
phase space of system (18), because in the extended symplectic space R^^ the 
phase space has measure zero. 

2.2. Existence and uniqueness of the motion. The Cauchy problem for the 
transport equation related to ODE (18) is as follows: 

(27) ^ = (^Ax - Bsign(^B,^{x)^ , v{x,0)=v{x). 

Our main result claims that in the phase space R^^ of the system (18) we can define 
a semiflow which is continuous, uniquely defined everywhere, and it is related to 
the transport equation (27) in a way the flow from Theorem 2 is related to (20): 
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Theorem 3. There exists a unique continuous semiflow (/'t(x), t > 0 such that 
v{x,t) — is the unique renormalized solution of the Cauchy problem for 

the transport equation (27). Every curve x{f) = (ft{x) is absolutely continuous and 
differential inclusion (18) holds. 

The proof is rather lengthy and can be found in Ref. [4]. 

This continuity implies, in particular, that the flow is defined uniquely every¬ 
where, although the control u{x) is defined uniquely only outside the hypersurface 
{{B,dp/dx)=0}. A similar phenomenon was discovered by LA. Bogaevskii [15] for 
the gradient differential equations x = —dffdx, where / is a nonsmooth convex 
function. 

The singular part of the right-hand side of (18) has the form —a{x)^, where 
a is a smooth nonnegative symmetric matrix, while f is a (nonsmooth) convex 
function. The quadratic form 

(28) (a(x)C,0 + 

is strictly positive, and the singular part of (18) is invariant under scaling x i—>■ 
Xx of the phase space. These facts allows us obtain differential inequalities for 
(a(a;)|j, |j) and {x, |y)^, where u is a solution of (27), and a priori bounds for the 
Lipschitz constant of v in any domain of the form {|(()t(a:)| > c}. More precisely, we 
need differential inequalities for a solution of regularized equation (27), where sign 
is replaced by a smooth function approximating sign in Li. In this way, we obtain 
the Lipschitz constants which do not depend on the accuracy of the regularization. 

Conclusion 

The generalized dry-friction control for a system of oscillators is especially inter¬ 
esting because of its asymptotic optimality, previously studied in Refs. [4-6]. As it is 
usual in the optimal control theory, this control is not defined uniquely everywhere 
and is a discontinuous function of the phase state of the system. Existence of the 
motion under control follows from the Filippov theory. However this theory does 
not answer the question of uniqueness of trajectories. It is shown that the existence 
and uniqueness issue of the motion under the generalized dry friction control can be 
resolved in the framework of DiPerna-Lions theory. It is important and interesting 
to explore similar issues for the optimal control. 

In real life, usage of a discontinuous control is accompanied with difficulties of 
various kinds. In the considered case, the discontinuous control in the form of a 
generalized dry friction arise naturally from the study of the asymptotic behavior 
of reachable sets of general linear systems. As noted above, in the proof of the main 
Theorem 3 we use an approximation of a singular law of motion by smoothing 
the function x i—>■ sign(a;). This quite common method, of course, can be useful in 
practical applications. 

We emphasize that the control law in the generalized dry friction form does 
not completely solve the problem of damping the oscillations of a system of linear 
oscillators. The method proposed in author’s papers [4-6] includes a combination of 
three control strategies: for large, intermediate and small values of the energy of the 
system. In the case of high energy, generalized dry friction is used. For intermediate 
values of the energy we use a similar control with a smaller amplitude. By decreasing 
the amplitude of the control we are reducing the size of the standstill zones. In a 
close vicinity of the equilibrium, this allows us to switch to a completely different 
type of control based on common Lyapunov functions. 

Thus, the result of the paper on the existence and uniqueness of the motion under 
the generalized dry friction is relevant at the first two stages of the motion to the 
terminal set, while at the last stage the issue has a different nature because of the 
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substantially different control law. From the viewpoint of asymptotic optimality, 
the motion at the first stage plays a decisive role. 
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